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EVALUATION  OF  PRECISION  AND  ACCURACY  IN  THE  CALIBRATION  OF  HYDROPHONES 


INTRODUCTION 

When  a  standardizing  institution  such  as  the  National  Bureau  of 
Standards  issues  a  certificate  for  a  standard,  it  attributes  to  that 
standard  a  numerical  value  together  with  a  statement  about  the  uncertainty 
of  this  value  [l] .  An  uncertainty  statement  is  necessary  in  calibration 
because  from  a  philosophical  viewpoint  "absolute  certainty  is  a  privilege 
of  uneducated  minds  and  fanatics"  [2].  The  standardizing  laboratory  never 
can  be  absolutely  certain  that  the  certified  value  falls  within  the  limits 
stated  in  the  certificate;  however,  from  a  practical  viewpoint,  it  is 
quite  justified  in  having  confidence  in  a  calibration  value  based  on  its 
knowledge  of  the  calibration  process  and  the  physical  nature  of  the  device 
calibrated.  This  confidence  generally  is  transmitted  in  the  form  of  an 
uncertainty  statement. 

The  primary  function  of  the  Underwater  Sound  Reference  Division  (USRD) 
of  the  Naval  Research  Laboratory  (NRL)  is  to  serve  as  the  Navy's  standard¬ 
izing  institution  for  underwater  electroacoustic  measurements.  Because 
stable  electroacoustic  devices  are  required  if  calibration  is  to  remain 
valid  for  a  reasonable  length  of  time  and  under  a  range  of  environmental 
conditions,  and  because  a  thorough  physical  knowledge  of  the  standards  Is 
necessary  to  good  calibration,  an  important  role  in  fulfilling  the  pri 
mary  mission  is  played  by  the  Standards  Branch  of  the  USRD  in  the  develop¬ 
ment  and  manufacture  of  stable,  transportable,  standard  underwater 
electroacoustic  transducers  for  Navy  activities  and  contractors. 

At  the  USRD,  knowledge  of  the  characteristics  of  underwater  electro¬ 
acoustic  transducers  and  of  the  calibration  process  has  led  to  the  belief 
that  the  common  calibration  procedures  employed  have  negligible  bias  and 
that  the  accuracy  associated  with  these  measurements  is  ±1  dB.  No  confi¬ 
dence  limits  are  specified,  and  the  accuracy  is  an  educated  guess  based 
on  experience.  The  purposes  of  this  study  are:  (1)  to  establish  valid 
confidence  limits  on  the  calibration  outputs  of  the  existing  calibration 
facilities;  (2)  where  possible,  separate  and  identify  systematic  and 
random  errors  associated  with  each  facility;  and  (3)  where  appropriate, 
recommend  changes  in  calibration  procedures,  or  equipment,  or  both,  to 
decrease  the  confidence  interval  associated  with  the  calibration 
measurements . 
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Because  the  accumulation  of  enough  data  to  establish  valid  confidence 
limits  is  a  time-consuming  process,  and  the  USRD  calibration  systems’ 
workloads  are  such  that  time  is  available  to  acquire  the  needed  data  only 
infrequently,  the  results  in  this  report  are  based  on  a  very  limited 
amount  of  data;  they  should  be  viewed  as  only  an  initial  look  at  a  pr:b- 
lem  that  will  require  a  more  thorough  system-by-system  analysis.  Results 
of  this  study  do  indicate  which  systems  require  the  most  immediate  atten¬ 
tion  and  which  can  continue  to  be  used  with  a  reasonable  degree  of 
certainty. 

As  measures  of  experimental  error,  precision  and  accuracy  as 
defined  by  Mandel  [l]  are  used.  Precision  is  defined  in  terms  of  its 
opposite,  imprecision.  Imprecision  is  the  amount  of  scatter  exhibited  by 
results  obtained  from  repeated  experimental  measurements.  There  will 
always  be  some  imprecision  associated  with  experimental  measurements 
because  of  lack  of  fineness  of  scale  on  some  measuring  devices,  although 
the  results  may  seem  to  be  precise.  The  accuracy  of  an  experimental 
measurement  can  be  lefined  in  terms  of  the  bias  or  systematic  error  of 
the  experimental  measurements  or  as  the  difference  between  the  true  value 
and  the  expected  value,  which  can  only  be  estimated  from  repeated  measure¬ 
ments.  From  these  definitions,  one  can  see  that  a  measurement  procedure 
can  be  considered  precise  but  not  accurate  or  accurate  but  not  precise. 
Using  these  definitions  avoids  some  of  the  philosophical  problems  asso¬ 
ciated  with  another  common  definition  of  accuracy  (that  is  concerned  with 
both  bias  and  precision)  and  is  admittedly  taking  the  easy  way  out  [3]. 

In  brief,  we  avoid  the  situation  in  which  a  procedure  having  a  small  bias 
and  a  high  precision  is  called  more  accurate  than  an  unbiased  procedure  of 
low  precision.  Uncoupling  these  terms  allows  one  to  distinguish  between 
procedures  in  which  repetition  can  reduce  the  confidence  interval  about 
an  unbiased  mean  and  those  in  which  repetition  beyond  a  certain  limit 
does  no  good  because  the  mean  is  biased. 

The  first  problem  one  encounters  in  establishing  confidence  intervals 
for  the  calibration  of  hydrophones  is  that  of  arriving  at  a  philosophi¬ 
cally  acceptable  definition  of  the  true  value  of  hydrophone  sensitivity  in 
terms  of  ideal  deterministic  parameters,  because  the  true  value  must  be 
known  before  one  can  say  how  much  error  is  involved  in  a  single  measure¬ 
ment.  This  problem  has  been  solved  to  some  extent  by  the  application  of 
electroacoustic  two-port  theory  in  reciprocity  calibration  [4].  Under 
ideal  conditions,  reciprocity  calibration  is  supposed  to  be  an  absolute 
method;  however,  several  conditions  must  be  satisfied  to  justify  the  con¬ 
clusion  that  a  free-field  reciprocity  calibration  results  in  an  unbiased 
or  absolute  value  of  hydrophone  sensitivity.  These  conditions  are  the 
existence  of  a  free  sound  field,  sufficient  space  to  permit  measurements 
in  the  far  field  of  the  source,  linearity,  and  infinite  signal-to-noise 
ratio  (SNR),  none  of  which  can  be  achieved  in  practice,  so  some  bias  or 
inaccuracy  in  the  measurement  is  inevitable.  This  paper  discusses  the 
tiases  to  be  expected  in  limited-water  reciprocity  calibration  and 
attempts  to  estimate  their  magnitude  or  minimize  their  effects.  Hydro¬ 
phones  calibrated  repeatedly  after  minimizing  biases  have  been  used  in 
this  study  as  primary  standards  having  "known"  calibrations  to  provide  an 
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initial  evaluation  of  the  practices  and  systems  used  at  the  USRD  in  the 
established  calibration  facilities - 

Early  in  this  investigation  it  was  recognized  that  the  evaluation  of 
system  accuracy  and  precision  does  not  provide  a  confidence  interval  and 
bias  statement  for  every  type  of  hydrophone  that  might  be  calibrated  on 
that  particular  system;  however,  some  general  statements  about  extreme 
confidence  intervals  and  bias  may  be  made  for  certain  categories  of 
hydrophones  such  as  the  NRL-USRD  standards  operated  below  resonance  for  a 
particular  system.  This  region  for  these  standards  is  of  primary  interest 
to  the  NRL-USRD  in  fulfilling  its  primary  mission  of  serving  as  the  stand¬ 
ardizing  institution  for  underwater  electroacoustic  measurements.  The 
accuracy  and  precision,  which  in  our  other  measurements  is  of  secondary 
interest,  depend  on  the  instrument  upon  which  the  measurements  are  being 
performed,  and  often  is  so  dependent  upon  the  instrument  that  additional 
charges  sometimes  should  be  made  to  the  customer  if  very  good  evaluation 
of  the  measurements  is  desired.  The  present  investigation  is  confined  to 
the  problem  of  specifying  accuracy  and  precision  of  standard  hydrophones 
at  frequencies  below  resonance. 


THEORY 

To  evaluate  the  precision  and  accuracy  in  calibration  of  hydrophones 
by  the  reciprocity  method  requires  an  understanding  of  electroacoustic 
reciprocity  and  calibration  practice.  MacLean  [4]  developed  the  reci¬ 
procity  method  for  the  absolute  calibration  of  electroacoustic  transducers. 
His  development  was  confined  to  a  lumped-parameter  treatment,  which  is  not 
generally  valid  when  the  motion  of  the  electroacoustic  transducer  is  not 
the  same  for  projecting  as  it  is  for  receiving  sound  waves.  Foldy  ard 
Primakoff  [5]  proved  that  the  reciprocity  method  was  independent  of  one 
particular  characteristics  of  the  transducer  so  that  the  differences  in 
motion  in  projecting  and  receiving  do  not  matter.  Their  rather  general 
proof  was  based  on  generalized  impedances.  Green's  functions,  and  both 
acoustic  and  electroacoustic  reciprocity.  Bobber  [6]  has  discussed  both 
the  electroacoustic  reciprocity  calibration  method  and  practice  thoroughly 
in  his  recent  book  on  underwater  electroacoustic  measurements.  Also,  an 
understanding  of  the  effect  of  resonances  and  radiation  impedance  on 
hydrophone  sensitivity  is  essential  in  evaluating  the  precision  and  accu¬ 
racy  of  calibration,  because  a  properly  designed  hydrophone  can  be  cali¬ 
brated  with  better  accuracy  than  can  a  poorly  designed  one. 

Here,  the  basic  definition  of  hydrophone  sensitivity  will  be  reviewed 
using  lumped-parameter,  electromechanical  two-port  theory  instead  of  the 
more  generalized  approach  of  Foldy  and  Primakoff.  Then  bias  in  the 
measurement  of  hydrophone  sensitivity  by  the  reciprocity  method  will  be 
discussed.  This  will  be  followed  by  a  review  of  the  statistical  concepts 
associated  with  the  acquisition  and  averaging  of  the  hydrophone  sensi¬ 
tivity  estimates. 


3 


Electromechanical  Two-Port  Theory  for  a  Hydrophone 

If  a  transducer  is  linear  and  reciprocal,  it  can  be  represented  as  a 
simple  electromechanical  two-port  device  as  illustrated  in  Fig.  1.  The 


Fig.  1.  Representation  of  an  electromechanical  two- 
port  device. 


two-port  equations  for  this  device  (assuming  the  vibratory  mode  is  the 
same  whether  projecting  or  receiving)  are  given  by 
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where  p  is  the  average  applied  pressure  over  the  area  A  of  the  mechanical 

port,  vn  is  the  average  velocity  normal  to  A,  i  is  the  electrical  current 

in  the  electromechanical  two  port,  e  is  the  voltage  across  the  electrical 
port,  and  the  impedance  parameters  are  defined  by  the  equations 
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The  quotient  of  open-ci  L-cuit  voltage  divided  by  the  force  on  area  A  is 
obtained  from  Eqs.  (1)  and  (2)  as 
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To  obtain  an  output  voltage  from  an  electromechanical  device,  a  force 
must  be  delivered  to  the  mechanical  port.  Perce  is  delivered  to  a  hydro¬ 
phone  through  vibration  of  the  medium.  The  mechanical  port  looks  out  upon 
the  medium  and  its  boundaries  as  a  mechanical  radiation  impedance  Z 

mr 

The  coupling  of  the  force  from  the  mechanical  radiation  impedance  to  the 

mechanical  port  probably  can  be  best  understood  in  terms  of  a  purely 

mechanical  system.  From  the  dual  mechanical  port  of  Fig.  1,  one  obtains 

Fig.  2;  force  Ap  is  transmitted  by  the  rod-connected  system  from  Z  to 

mr 

Z^  when  the  electrical  port  is  open-circuited.  Note  that  the  rules  appli¬ 
cable  to  a  mechanical  system  are  not  the  same  as  those  for  an  electrical 
circuit.  Impedances  in  parallel  add  because  they  are  defined  as  the 
through  variable  (force)  divided  by  the  across  variable  (velocity) .  Using 
the  mechanical  analog  of  Norton's  equivalent  circuit,  one  obtains  the 
mechanical  system  shown  in  Fig.  3,  where  F  is  the  average  force  that 


Fig.  2.  Mechanical  system  for 
an  open-circuited  electro¬ 
mechanical  transducer. 


Fig.  3.  Equivalent  mechanical 
system  using  the  analog  of 
Norton's  equivalent  circuit. 


arises  when  the  diaphragm  of  the  transducer  is  blocked.  The  blocked  force 
is  equal  to  the  product  of  the  area  and  the  average  blocked  pressure  p^. 
From  Fig.  3  it  is  deduced  that 
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The  average  blocked  pressure  is  related  to  the  free-field  pressure  pf  by 
the  equation 


(9) 
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where  D  is  the  diffraction  constant  of  the  transducer.  The  receiving 
sensitivity  of  the  transducer  is  obtained  from  Eqs.  (7)/  (8),  and  (9): 


Bias  in  Reciprocity  Measurement  of  Hydrophone  Sensitivity 

The  reciprocity  method  calibration  equation  for  free-field  voltage 
sensitivity  given  by  Bobber  [6]  in  terms  of  deterministic  quantities  or 
true  values  is 

P  A 

TH  PH 

M  =  -  • 

ePTiT 

where  P,  T,  and  H  are  symbols  for  the  projector,  transducer,  and  hydro¬ 
phone,  respectively;  the  e's  are  open-circuit  voltages  under  noise-  and 
interference-free  conditions,  and  the  first  subscript  refers  to  the  sound 

emitter  while  the  second  refers  to  the  receiver;  i  is  the  transducer 

T 

current;  d  is  the  measurement  distance;  p  is  the  water  density;  and  f  is 
the  frequency.  Equation  (11)  gives  the  sensitivity  in  volts  per  micro¬ 
pascal  if  the  e's  are  in  volts,  i^  in  amperes,  d  in  centimeters,  p  in  grams 

per  cubic  centimeter,  and  f  in  hertz.  In  measuring  the  quantities  em„, 

e„  .  and  e  ,  the  distance  d  is  measured  each  time  so  that  each  measured  e 
PH  PT 

has  a  voltage  measurement  error  and  a  distance  measurement  error  asso¬ 
ciated  with  it.  In  free-field,  far-field  measurements,  the  functional 
dependence  of  the  true  values  of  the  open-circuit  voltages  is  given  by 

e  =  eg/d,  (12) 

where  eg  is  the  voltage  that  would  appear  at  distance  dg-  Because  measure¬ 
ments  are  made  in  the  presence  of  noise,  an  estimate  of  e  for  distance  d 
is  given  by 

e  -  (e2  +  cn2)Sj(d/d) 

where  the  circumflex  represents  an  estimated  quantity  for  a  random  varia¬ 
ble,  and  an  is  the  root-mean-square  (rms)  value  of  an  assumed  zero-mean, 

stationary,  Gaussian  noise,  and  e  is  an  rms  value.  Here  it  was  assumed 
also  that  signal-to-noise  ratio  does  not  change  with  small  changes  in 
distance.  An  estimate  of  M  can  now  be  written  as 
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In  Eq.  (14)  all  of  the  right-hand  quantities  with  the  circumflex  can  be 
considered  independent.  The  question  of  distribution  and  bias  in  the 
estima . .  n  of  M  now  arises .  Provided  that  the  errors  are  small  in  compari¬ 
son  with  the  quantities  being  measured,  the  approximate  error  in  computing 
the  sensitivity  can  be  found  by  taking  the  variation  of  Eq.  (14) .  Doing 
this,  dividing  the  result  by  Eq.  (14),  and  replacing  the  estimated  quanti- 
tit.ies  in  Eq.  (14)  by  their  deterministic  counterparts,  one  obtains 
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for  the  relative  uncertainty  in  the  derived  measurement  of  hydrophone 
sensitivity  by  the  reciprocity  method.  For  a  measurement  process,  rela¬ 
tive  uncertainty  is  a  random  variable.  To  compute  a  confidence  interval 
for  the  measurement  process,  it  is  necessary  to  know  or  to  approximate  the 
probability  density  function  for  the  relative  uncertainty  of  the  process. 
An  uncertainty  is  associated  also  with  the  measurement  of  the  probability 
density  function  for  the  relative  uncertainty.  For  these  reasons,  rela¬ 
tive  uncertainty  often  is  considered  v>  be  a  zero-mean,  Gaussian  variable, 
but  without  adequate  justification.  For  the  case  of  free-field  measure¬ 
ments  made  with  freely  suspended  electroacoustic  transducers,  one  can 
argue  that  the  assumption  is  valid  under  certain  conditions.  For  this 
case,  the  only  factors  contributing  Significantly  to  the  relative 
uncertainty  of  hydrophone  sensitivity  are  the  relative  uncertainties  of 
voltage  and  distance  values.  The  voltage  estimates  can  be  assumed  to 
arise  from  a  Gaussian  process  if  noise  is  stationary  and  if  the  time  con¬ 
stant  of  the  measuring  instrument  is  large  in  comparison  with  the  recip¬ 
rocal  of  the  bandwidth  of  the  interfering  noise  because  of  the  central 
limit  theorem.  Distance  can  be  measured  between  the  centers  of  the  rig¬ 
gings  for  the  projectors  and  hydrophones,  but  the  distance  between  the 
acoustic  centers  of  these  instruments  may  vary  randomly  and  may  also  con¬ 
tain  a  bias  distance.  If  these  instruments  are  cylindrical  and  can  be 
suspended  by  their  cables,  one  may  argue  that  the  distance  can  be  meas¬ 
ured  at  the  air-water  interface  and  that  distance  bias  can  be  ignored. 
Repeated  rigging  then  probably  would  result  in  the  relative  uncertainty 
of  distance  being  an  unbiased  Gaussian  variable.  Because  the  sum  or  dif¬ 
ference  of  Gaussian  variables  results  in  a  Gaussian  variable,  the  rela¬ 
tive  uncertainty  of  hydrophone  sensitivity  as  well  as  the  estimates  of 
hydrophone  sensitivity  can  be  considered  a  biased  Gaussian  variable. 
Hydrophone  sensitivity  is  commonly  expressed  on  a  logarithmic  basis  refer¬ 
enced  to  1  pbar  or  to  1  pPa.  It  is  well  known  that  the  average  of  the 
logs  always  will  be  less  than  or  equal  to  the  log  of  the  average.  Hershey 
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[7]  has  shown  that  when  the  ratio  of  the  standard  deviation  to  the  mean 
a/y  is  less  than  1/10  (fluctuation  -  ±1  dB) ,  then  there  is  less  than  a 
0.05-dB  difference  between  the  average  of  the  logs  and  the  log  of  the 
average  for  a  Gaussian  process  when  the  average  is  taken  over  fewer  than 
30  values. 

Careful  attention  to  the  output  voltage  of  transducer  or  hydrophone 
will  ensure  that  the  readings  are  as  free  of  bias  as  possible.  There  are 
four  apparent  sources  of  bias  in  the  voltage  readings:  (1)  ambient  noise, 
(2)  crosstalk,  (3)  multipath  transmission,  and  (4)  bias  in  reading  and  in 
voltmeter  calibration.  The  effect  of  ambient  noise  is  as  indicated  in 
Eq.  (13) .  Keeping  the  signal-to-noise  ratio  above  20  dB  reduces  noise 
bias  to  less  than  0.05  dB  on  the  average.  Crosstalk  can  be  a  more  serious 
problem  and  should  be  kept  40  dB  below  signal  level  to  limit  bias  to  less 
than  0.05  dB.  Bias  in  voltmeter  calibration  and  reading  is  virtually 
eliminated  by  choice  of  a  good  voltmeter.  Effects  of  multipath  trans¬ 
mission  caused  by  lack  of  true  free-field  conditions  must  be  taken  into 
account . 

When  the  free-field  sensitivity  of  a  hydrophone  is  measured  in  a 
bounded  environment  where  multipath  transmission  occurs,  a  change  of 
sensitivity  arises  from  the  change  in  mechanical  radiation  impedance  and 
in  bias  caused  by  multipath  acoustic  pressures.  As  an  example  of  how  a 
boundary  affects  radiation  impedance,  consider  a  small  spherical  trans¬ 
ducer  P  of  ^dius  a  operated  near  an  air-water  interface  as  illustrated 


Fig.  4.  Geometry  for  projector  P,  which  estab¬ 
lishes  pressure  Po  in  the  vicinity  of  an  air- 
water  interface. 


in  Fig.  4.  Assume  h  »  a  and  the  wavelength  \  radiated  is  large  .in  com¬ 
parison  with  a  (monopole  source) .  Let  the  spherical  radiator  pulsate 
with  radial  velocity 


v 

r 


(16) 


where  w  is  the  angular  frequency,  and  vrQ  is  the  radial  velocity  at  the 
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surface  of  the  sphere.  When  the  amplitude  of  the  vibrations  is  small  in 
comparison  with  a,  a  «  d.  Under  free-field  conditions,  the  acoustic 
pressure  at  point  Q,  distance  d  from  the  center  of  the  sphere,  is  given  by 


p  (d/t)  =  -  k(d-a)],  (17) 

d 

where  pg  is  the  pressure  amplitude  at  the  surface  of  the  sphere  and  k  is 
the  wave  number.  The  specific  acoustic  radiation  impedance  is  obtained 


(ZSr)  ff  =  pff(a't)/vr  =  VVrO  (18) 

for  the  free-field  case.  Because  the  wavelength  is  assumed  to  be  large  in 
comparison  with  the  radius  of  the  sphere,  the  spacific  acoustic  radiation 
impedance  can  be  shown  to  be 


pck2a2  +  jpcka 


(Zsr*ff 


1  +  k2a2 


where  c  is  the  speed  of  sound  in  water  and  p  is  the  mean  density.  Now 
consider  the  same  spherical  radiator  pulsating  with  the  same  radial 
velocity  but  located  at  depth  h  below  the  air-water  interface.  The  acous¬ 
tic  pressure  at  a  point  d  distant  from  the  center  of  the  sphere  and  at  the 
same  depth  is  given  by 


P,..(d,t)  =  pa 

aw  s 


^-jkCd-a)  -jkCdj-a)]  . 


where  the  subscript  aw  refers  to  the  air-water  interface,  and 


dj  =  (4h2  +  d2)*5. 


The  specific  acoustic  radiation  impedance  for  this  case  is  given  by 


Paw(a,t)  e_2jkh  (22 

(Z  )  2  — -  =  (Z  )  1 - . 

sr  aw  v  sr  ff  2kh 

r 

Specific  acoustic  impedances  are  converted  to  mechanical  impedances  by 
multiplying  by  the  area  of  the  sphere.  The  mechanical  radiation  imped¬ 
ances  are  approximately 


^Zmr^  ff 


4TTa2(puka2  +  jwpa) 
1  +  k2a2 
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(24) 


The  effect  of  the  mechanical  radiation  impedance  in  the  presence  of  an  air- 
water  interface  on  hydrophone  sensitivity  is  evident  in  the  universal  elec¬ 
trical  equivalent  circuit  of  a  Class  I  transducer  (Fig.  5),  which  is  given 


Fig.  5.  Electrical  equivalent  cir¬ 
cuit  for  open-circuited  Class  I 
transducer. 


by  Fischer  [8]  (some  notational  changes  have  been  made) .  With  the  circuit 
viewed  as  a  voltage  divider, 
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■ff 


4>AD 


1  +  j(42C  (Z. 
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+  Z  ) 
mr 


(25) 


is  obtained  for  hydrophone  sensitivity,  where  <J>  is  the  electromechanical 
conversion  factor.  The  effect  of  radiation  impedance  on  sensitivity  is 
negligible  far  below  resonance,  where  sensitivity  is  controlled  primarily 
by  the  mechanical  compliance  C^  of  the  transducer  element  and  is  essen¬ 
tially  constant  with  frequency  provided  that  D  is  constant  with  frequency 
over  this  region.  Sensitivity  is  most  seriously  affected  by  interfering 
boundaries  at  resonance  because  the  radiation  mass  can  become  comparable 
with  the  vibrating  mass  of  the  transducer  element  and  radiation  reactance 
is  zero.  For  example,  the  vibrating  mass  of  a  spherical,  rigidly  backed 
ceramic  shell  vibrating  in  the  pulsating  mode  is  about  half  that  of  the 
shell,  while  the  radiation  mass  is  three  times  that  of  the  displaced 
water  (for  the  free-field  case).  For  other  hydrophone  configurations,  the 
effect  of  water  mass  loading  on  sensitivity  decreases  with  increasing 
directivity. 


An  example  of  another  source  of  bias  in  the  calibration  of  hydrophones 
by  reciprocity  is  that  produced  by  the  interference  pressure  at  an  air- 
water  interface  (excluding  all  other  effects) .  The  signal  reflected  from 
the  interface  and  the  direct  signal  sum  constructively  and  destructively 
to  produce  a  cyclic  amplitude  as  a  function  of  frequency,  which  is  known 
as  the  interference  cycle  because  it  is  periodic  with  frequency.  When  the 
estimates  of  hydrophone  sensitivity  are  made  over  an  interference  cycle, 
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the  bias  B(M)  relative  to  free-field  calibration  is 


B  (M)  M 
_ aw 

M  M 


(26) 


where  M  denotes  the  sensitivity  measured  near  the  air-water  interface, 
aw 

The  reciprocity  equation  gives 


B(M) 


M 


(<Waw  (eTHJaw  ®PT 


e  e  (e  ) 

PH  TH  PT  aw 


-  1. 


(27) 


If  it  is  assumed  that  far  below  hydrophone  resonance  the  sensitivities  for 
free-field  and  near  the  air-water  interface  are  the  same,  then 


(£Waw  "  M(Pu) 


H'aw' 


(28) 


and 


epH  =  MpH'  (29) 

etc . ,  where  (p  )  and  p  are  the  pressures  at  the  hydrophone  near  an  air- 
H  aw  n 

water  interface  and  in  a  free  field,  respectively.  Application  of  these 
relationships  to  Eq.  (26)  gives 


B(M) 
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(<Waw 


'TH 
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(30) 


The  actual  bias  obtained  will  depend  upon  the  detection  method  used  in 

obtaining  e_It  and  (e_„) 

3  TH  TH  aw 

For  rms  detection  in  a  free  field,  let  the  instantaneous  pressure 
field  at  the  hydrophone  due  to  the  reciprocal  transducer  be 


p  (d,t)  «  ^n(a3t-~.  k<3) 
yH  '  d 


(31) 


After  squaring,  time  averaging,  and  taking  the  square  root  of  this  expres¬ 
sion,  one  obtains 

PH  «  l/d^2.  (32) 

In  the  presence  of  an  air-water  interface,  the  pressure  at  the  hydrophone 
is 
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sin (tot  -  kd)  sin[tot  -  k(4h2  +  d2)  5] 
d  (4h2  +  d2)** 


Squaring,  time  averaging  over  an  integral  number  of  periods,  and  taking 
the  square  root  of  this  equation  gives 


’<*  [l 

[p  (d,t)]2  dt  «  - + - 

H  aw 


cos  k[(4h2  +  d2)*5  -  d]]*5 


2d2  2  (4h2  +  d2)  d(4h2  +  d2)*5 


.  (34) 


Averaging  Eq.  (34)  over  an  interference  cycle  produces 

r2Ttr 


8h2  +  4d2 


cos  kAx 


(Pli)aw  2ir  Q  [4d2(4h2  +  d2)  d(4h2  +  d2)^ 

//iv.2  j.  aIs**  _  j 


d(kAx)  , 


where  Ax  =  (4h2  +  dz)  -  d. 

The  ratio  of  the  rms-detected  pressures  averaged  over  an  interference 
cycle  is  then 


<Vaw  1  K  +  2dT 


2ir  4h2  +  d2 


2d(4h2  +  d2)^  ^ 

1 - cos  kAx  d(kAx).  (35) 

4h2  +  2d2 


Let  b  =  d/h.  Then  Eq.  (36)  becomes 

(  f2lT  r 


<Vaw  1  4 


+  2b2  ** 


2tt  14  +  b^ 


2b  (4  +  b2)*2  55 

X - cos  kAx  d  (kAx) . 

4  +  2b2 


Expanding  the  integral  into  a  binomial  series  and  integrating,  one  obtains 
(pH}aw 


4  +  2b2  «  [: 

- —  .  1  _  l  - 

4  +  b2  n=l 

. 

for  all  values  of  b. 


Cl* 5  •••  (4n  -  3)][l*3  •••  (2n  -  1)] 

*  ' 

b2  (4  +  b2) 

(2n) !  [2*4  •••  (2n) ] 

(4  +  2b2)2 

. 

Equations  (30)  and  (38)  have  been  used  to  determine  the  bias  in  cali¬ 
brating  a  small  (a  «  X,  h  >>  a)  nonresonant  hydrophone  by  the  reciprocity 
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method  in  the  vicinity  of  an  air-water  interface.  Table  I  gives  the 
results  of  this  calculation  for  several  separation-to-depth  ratios.  For 
most  practical  purposes  a  separation-to-depth  ratio  of  1/4  can  be  taken  to 
give  an  unbiased  calibration  when  the  bottom  surface  is  remote. 

Table  I.  Bias  in  reciprocity  cali¬ 
bration  near  an  air-water  interface 
with-  c-w  projected  energy;  b  =  d/h; 

(Pr) aw  is  pressure  at  the  hydrophone 

near  the  interface;  p  is  free-field 

H 

pressure  at  the  hydrophone;  B(M)  is 
the  relative  bias;  M  is  the  free- 
field  sensitivity;  d  is  the  distance 
from  the  center  of  a  spherical  radia¬ 
tor  to  the  measuring  point;  and  h  is 
the  distance  from  the  radiator  to 
the  interface. 


b 

<Vaw 

PH 

B  (M) 

M 

B(M) 

(dB) 

3.00 

1.190 

0.091 

+0.76 

2.00 

1.133 

0.064 

+0.54 

1.50 

1.0b'.' 

0.044 

+0.37 

1.00 

1.051 

0.025 

+0.21 

0.50 

1.028 

0.014 

+0.12 

0.25 

1.004 

0.002 

+0.02 

For  peak  detection  under  free-field  conditions  with  negligible  noise, 

PH  «  1/d.  (39) 

In  the  presence  of  an  air-water  interface,  Eq.  (33)  applies  for  the  pressure 

field  at  the  hydrophone.  Let  dj  =  (4h2  +  d2)*2.  Then  Eq.  (33)  can  be 
written  as 

(di  coskd -d  coskdi)  sin  lot  +  (d  sinkdi  -di  sinkd^  cosiot  , 

[p„(a,t)]...  -  — - - - - - .  (40> 


aw 


ddi 


The  amplitude  of  the  wave  represented  by  Eq.  (40)  is  the  square  root  of 
the  sum  of  the  squares  of  the  quadrature  components,  or 
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(41) 


[d|  +  d2  -  2ddj  cos  k(dj  -  d) ] * 


Averaging  Eg.  (si)  over  an  interference  cycle  produces 


(p  )  K 

fH  aw 


f2ir  j, 

1  [d?  +  d2  -  2dd!  cos  kAx]1 


d(kAx)  . 


The  ratio  of  the  peak-detected  pressures  averaged  over  an  interference 
cycle  is  then 


Vaw  1  .  4  -  2b2 

P..  2tt  14  +  b2 


,  i.  )Jj 

2b  (4  +  b2) ^ 

1  -  -  cos  kAx  d(kAx).  (43) 

4  +  2b2 


Equations  (37)  and  (43)  are  identical,  which  shows  that  the  biases  in 
reciprocity  calibration  for  rms-  and  peak-detection  near  an  air-water 
interface  are  the  same  when  there  is  no  noise.  This  result  is  as 
expected . 


Statistical  Concepts  for  the  Analysis  of  Hydrophone  Sensitivity 

The  proposition  has  been  put  forth  that,  under  certain  well  specified 
conditions,  reciprocity  calibration  leads  to  a  set  of  estimates  of  hydro¬ 
phone  sensitivity  that  are  unbiased  and  that  are  Gaussianly  distributed. 
Here  a  review  of  the  statistical  parameters  needed  to  specify  the  confi¬ 
dence  coefficient  associated  with  measurements  from  a  Gaussian  process 
will  be  given. 

Associated  with  an  infinite  Gaussian  population  of  estimates  are  two 
parameters  necessary  to  specify  the  Gaussian  probability  density  function 
for  that  population.  These  are  the  mean  value  M  and  the  variance  a^2 
defined  by  the  equations 

!  N 

M  =  lim  -  l  M  (44) 

N-*»  N  i=l 


1  N 

a2  =  lim  -  l  (M.  -  M)2. 

**  *T  *  _ *1  ^ 


N-*»  N  i=l 


The  Gaussian  probability  density  function  is  given  by 
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P  (M)  = 


1 


(46) 


p(M)  =  -  exp[-(M  -  M)2/2cr2J. 

It  is  well  known  that/  when  dealing  with  a  finite  number  of  estimates,  the 
average  of  the  estimates  is  a  better  estimate  in  a  probabilistic  sense 
than  is  a  single  estimate.  The  average  value  from  a  finite  set  of  esti¬ 
mates  taken  from  a  potentially  infinite  set  is  defined  by  uhe  equation 


1  N 

M  =  -  [  M  , 


N  i=l 


and  an  unbiased  estimate  of  the  standard  deviation  is  defined  by  the 
equation , 


5  =  - -  l  (M  -  5)2  . 

N  -  1  i=l  1 

The  estimate  for  the  average  of  N  values  of  M  is  itself  an  unbiased 
Gaussian  variable  with  a  standard  deviation  (standard  error) 


5m  =  VN  •  (49) 

The  confidence  interval  about  the  mean  associated  with  a  single  estimate 
is  written  as 

(M  -  A0m)  <  M  <  (M  +  A  a)  ,  (50) 

where  A  is  a  constant  to  be  determined.  The  confidence  interval  for  the 
average  of  the  estimates  is  written  as 

(M  -  cS./fT*)  <  M  <  (M  +  CO  /N*5)  ,  (51> 

where  C  is  a  constant  to  be  determined.  The  choice  of  the  constants  A 
and  B  depends  upon  the  value  of  the  confidence, coefficient  desired.  There 
is  a  difficulty  with  the  interval  of  Eq.  (50)  in  that  cr  is  not  known  for 

M 

measurement  of  hydrophone  sensitivity  but  can  only  be  estimated  by  repeated 

A  _  *  H 

measurements  as  cr Both  M  and  a  /N  in  Eq.  (51)  are  random  variables. 

M 


From  Eq.  (51) ,  one  defines 


M  -  M 
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The  random  variable  defined  by  the  quantity  inside  of  the  absolute  value 
signs  is  a  random  variable  which,  for  a  Gaussian  M,  follows  Student's  t 
probability  density  functions,  a  family  o*  functions  that  depend  upon  N. 
Like  the  Gaussian  distribution,  the  Student  t-distribution  is  symmetrical 
about  the  mean.  As  N  becomes  large,  the  Student  t-distribution  approaches 
the  Gaussian  distribution.  The  Student  t-distribution  can  be  used  to  com¬ 
pute  the  confidence  interval  associated  with  the  average  of  the  estimates 
and  has  been  tabulated  for  various  values  of  N  and  confidence  coefficients. 
The  confidence  coefficient  is  defined  as  the  area  under  the  Student's  t 
probability  density  function  covered  by  the  confidence  interval.  Because 
the  true  mean  value  is  fixed,  the  confidence  coefficient  expresses  the 
probability  that  the  random  confidence  interval  will  bracket  the  mean. 

What  should  one  do  if  the  estimates  of  hydrophone  sensitivity  do  not 
follow  a  Gaussian  distribution?  This  should  cause  no  real  alarm  (as  long 
as  the  bias  in  the  measurement  procedure  is  minimized)  because  of  the 
central  limit  theorem.  The  central  limit  theorem  is  expressed  as  follows: 
Given  a  population  of  values  with  a  finite  variance,  if  We  take  independ¬ 
ent  samples  from  this  population,  all  of  size  N,  then  the  population 
formed  by  the  averages  of  these  samples  will  tend  to  have  a  Gaussian  dis¬ 
tribution,  regardless  of  what  the  distribution  is  of  the  original  popula¬ 
tion;  the  larger  N,  the  greater  will  be  this  tendency  toward  the  Gaussian 
distribution  [l].  In  practice  N  does  not  have  to  be  very  large  for  the 
distribution  of  the  averages  to  be  nearly  Gaussian,  so  Eq.  (51)  will  apply 
relatively  well.  Also,  the  Student  t-distribution  can  be  used  to  express 
the  confidence  coefficient.  Because  there  are  some  differences  in  the 
Gaussian  distribution  and  the  distribution  of  the  averages,  high- 
confidence  coefficients  should  not  be  used  unless  N  is  fairly  large. 

Confidence  intervals  can  also  be  constructed  for  bias  estimation. 

The  bias  estimate  B  is  defined  by 

B  =  M  -  M.  (53) 

The  variance  of  §  is  the  same  as  the  variance  of  M,  since  M  is  fixed,  so 
the  confidence  interval  for  B  is 

(B  -  COg/W5)  <  B  <  (B  +  Cob/nS  ,  (54) 

where  5  is  the  average  of  the  B's.  The  interval,  like  that  for  the  confi¬ 
dence  interval  for  the  mean,  can  be  found  from  tables  of  Student's 
t-distribution . 


MEASUREMENTS  AND  ANALYSIS 

In  evaluating  the  accuracy  and  precision  of  existing  calibration 
facilities  at  the  USRD,  the  approach  chosen  was  to  calibrate  two  hydro¬ 
phones  by  reciprocity  at  frequencies  far  below  resonance  using  a  minimum 
of  equipment.  These  two  hydrophones  were  calibrated  with  high  accuracy 
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and  precision  several  times,  and  the  average  values  were  chosen  as  the 
"true"  values.  Then  the  same  hydrophones  were  calibrated  three  times  at 
each  calibration  facility  by  normal  procedures.  It  was  recognized  that 
three  calibrations  were  too  few  to  allow  good  evaluation,  but  workload 
requirements  prevent  making  a  large  number  of  measurements  at  each 
facility.  This  section  describes  the  measurements  and  gives  an  analysis 
of  them. 


Accurate  and  Precise  Calibration  of  Two  Standard  Hydrophones 

The  two  USRD  standard  hydrophones  used  in  this  study  were  the  F37 
serial  A40  and  the  F50  serial  21.  Drawings,  photographs,  and  typical 
sensitivity  curves  for  these  two  hydrophones  are  shown  in  the  Appendix. 
These  two  hydrophones  were  chosen  because  both  are  cylindrical  and  they 
can  be  suspended  by  their  cables  without  rigging  to  interfere  with  their 
calibration.  Also,  their  characteristics  are  known  to  be  relatively  con¬ 
stant  as  a  function  of  temperature  and  pressure.  At  10  kHz  and  lower 
frequencies,  both  hydrophones  can  be  considered  omnidirectional  in  the 
horizontal  plane  because  their  diameters  are  small  in  comparison  with 
a  wavelength. 

A  decision  was  made  to  attempt  to  calibrate  the  standard  hydrophones 
in  such  a  manner  as  to  ensure  that  a  confidence  interval  of  ±0.1  dB 
about  an  average  value  would  contain  the  true  mean  with  a  confidence 
coefficient  of  0.95  (95%  certain  to  contain  the  true  mean).  Figure  6  is 
a  block  diagram  of  the  system  used  in  the  reciprocity  calibration  of  the 
standard  hydrophones.  With  the  Fluke  voltmeter  and  the  Pearson  pulse  cur 
rent  transformer,  current  can  be  read  to  within  0.01  dB.  As  the  current 
transformer  rating  specifies  a  frequency  response  of  1  Hz  to  35  MHz  at 


Projector  Hydrophone 


Fig.  6.  Block  diagram  of  simple  calibration  system. 


the  -3  dB  points  and  an  output  voltage/ampere  of  0.1  with  a  range  of  0  to 
±0.1  dB  accuracy,  the  maximum  hydrophone  sensitivity  bias  expected  from 
the  current  measurement  bias  is  -0.05  dB,  which  can  be  deduced  from 
Eq.  (15) .  Measurement  of  open-circuit  hydrophone  voltages  to  within 
0.1  dB  is  possible  under  some  conditions.  The  environment  in  which  open- 
circuit  hydrophone  output  is  measured  determines  how  much  this  measured 
value  will  fluctuate.  In  Lake  Gem  Mary  at  the  USRD,  these  fluctuations 
were  so  serious  and  the  magnitude  of  the  boundary  interference  so  large 
that  an  unreasonably  large  number  of  calibrations  would  have  been 
required  to  obtain  the  desired  confidence  interval  with  the  simple  system 
shown  in  Figure  6.  Measurements  made  at  the  USRD's  Leesburg  Facility  with 
the  simple  system  and  with  the  transducers  at  about  mid-depth  revealed 
that  the  hydrophone  output  voltage  fluctuated  less  than  1%  under  most  con¬ 
ditions.  Exceptions  occurred  when  the  wind  blew  hard  and  when  fish  were 
in  the  vicinity  of  the  transducers.  Also,  wider  fluctuations  occurred 
when  the  signal-to-noise  ratio  was  insufficient  and,  at  the  lower  fre¬ 
quencies,  from  interference  by  60  Hz  and  its  harmonics.  Data  were  not 
recorded  when  fluctuations  exceeded  2%  of  the  voltmeter's  reading  and  then 
only  when  the  fluctuation  was  obviously  due  to  wind  conditions,  which 
caused  a  very  slow  fluctuation,  probably  because  of  a  lake  seiche.  One 
should  recognize  that  the  sound  pressure  field  is  changing  under  these 
conditions  even  though  no  measurable  change  in  projector  current  may  be 
observed,  because  the  radiation  impedance  may  be  small  in  comparison  with 
other  impedances  that  limit  the  current. 

The  two  hydrophones  were  calibrated  sequentially  in  a  single  reci¬ 
procity  calibration  procedure  with  both  the  low-cutoff  and  high-cutoff 
frequencies  of  the  Krohn-Hite  filter  set  on  the  frequency  at  which  the 
calibration  was  to  be  made.  Frequencies  chosen  were  0.500,  1.000,  2.000, 
5.000,  and  10.000  kHz.  Measurements  were  made  (a)  at  the  depth  21.06  m 
and  separation  2.00  m  with  two  F37's  as  the  projector  and  reciprocal 
transducer,  and  (b)  at  the  depth  17.36  m  and  separation  1.00  m  with  two 
J9's.  The  interference  cycles  for  the  air-water  interface  for  these  two 
set-ups  are  37.3  and  44.5  Hz,  respectively.  Twelve  frequencies  in  3-  and 
4-Hz  increments  about  the  chosen  frequencies  were  used  to  smooth  out  the 
interference  cycle  effects.  The  F37  and  J9  transducers  were  checked  for 
linearity  and  reciprocity.  They  were  found  to  be  both  linear  and  recip¬ 
rocal  to  within  ±0.04  dB  over  the  frequency  and  amplitude  ranges  in  which 
they  were  used.  A  distance  check  was  made  to  ascertain  how  accurately 
distance  was  being  maintained.  This  was  done  by  measuring  1-  and  2-m 
separations  on  two  free-hanging  F37's  and  averaging  the  received  hydro¬ 
phone  voltages  over  the  12  frequency  values  used  in  an  interference  cycle. 
The  results  of  these  measurements  show  that  the  distance  is  maintained  to 
within  ±1%  of  the  measured  distance  of  1  m,  which  makes  the  distance  error 
possibly  the  main  source  of  error  in  these  measurements.  The  error- is 
larger  with  the  F50  because  it  is  lighter  and  therefore  is  more  subject 
to  cable  deflections.  Nevertheless,  distance  is  essentially  unbiased 
when  it  .is  measured  repeatedly  and  averages  are  taken.  Independence  of 
the  measurements  was  assured  by  making  a  reciprocity  measurement  at  the 
12  frequency  points  comprising  one  interference  cycle.  All  measurements 
for  that  cycle  (that-  is,  three  voltages  and  the  current  for  each  point  in 


the  cycle)  were  made  at  the  same  instrument  settings  to  minimize  the 
effect  of  electronics  on  calibration. 


Figures  7  and  8  show  typical  plots  of  the  hydrophone  sensitivities 
calculated  by  using  the  free-field  reciprocity  parameter  with  J9's  used 
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fig.  7.  Unsmoothed  sensi¬ 
tivity  plots  for  F37  and  F50 
hydrophones . 


as  projector  and  reciprocal  transducer.  These  plots  show  the  effects  of 
boundary  interference.  It  was  calculated  that  the  interference  cycle  due 
to  the  air-water  interface  alone  was  about  44  Hz  and  that  the  peak-to-peak 
sensitive ty  variation  was  about  0.3  dB.  The  maximum  observed  variation 
resulting  from  the  interference  cycle  and  measurement  error  combined  for 
any  calibration  run  was  0.8  dB;  the  average  variation  was  about  0.5  dB. 
This  fact,  as  well  as  the  shapes  of  the  plots  in  Figures  7  and  8,  shows 
that  boundaries  other  than  the  air-water  interface  affected  the  measure¬ 
ments.  Averaging  over  the  12  points  taken  in  the  interference  cycle  is 
considered  to  provide  an  uabiased  estimate  of  hydrophone  sensitivity. 

Table  II  shows  the  results  of  sensitivity  measurements  made  on  hydrophones 
F37  serial  A40  and  F50  serial  21. 
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Table  II  shows  that,  with  15  independent  estimates  of  the  F37  hydro¬ 
phone  sensitivity  averaged  over  12  points  of  the  interference  cycle,  the 
maximum  spread  of  the  data  is  0.18  dB.  The  diffraction  constant  for  the 
F37  should  be  nearly  that  for  an  infinite  cylinder.  Using  the  equation 
for  the  diffraction  constant  giver,  by  Henriguez  [9],  the  effect  of  dif¬ 
fraction  at  10  kHz  is  estimated  to  reduce  the  sensitivity  by  about  0.1  dB 
from  that  at  lower  frequencies.  Taking  the  12  estimates  below  10  kHz  as 
being  estimates  for  the  same  true  value,  one  obtains  -204.71  dB 
re  1  V/uPa  for  their  average.  The  standard  deviation  of  these  points  is 
0.062  dB,  making  the  standard  error  of  the  average  5-  to  be  approximately 

0.02  dB.  From  Student's  t-distribution  table,  the  confidence  interval 
for  a  confidence  coefficient  of  0.95  is  M  i  0.04  dB.  This  means  that  the 
true  value  of  sensitivity  for  the  F37  serial  A40  lies  within  the  range 
-204.75  to  -204.67  dB  with  95%  certainty  if  the  measurement  is  unbiased 
as  hypothesized.  However,  some  small  positive  bias  is  sure  to  exist  so 
that  the  wider  range  -204.80  to  -204.67  dB  with  95%  certainty  is  not 
unreasonable. . 
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Table  II.  Simple  system  calibration  results;  estimate  of 
free-field  voltage  sensitivity  M  in  dB  re  1  V/pPa. 


Hydro-  Freq 

phone  (kHz)  J9  Reciprocity  F37  Reciprocity  M 


0.300 

-204.73 

-204.77 

-204.7 

1.000 

-204.73 

-204.67 

-204.7 

2.000 

-204.65 

-204.79 

-204.62 

-204.69 

-204.7 

5.000 

-2d4. 79 

-204.77 

-204.62 

-204.70 

-204.7 

10.000 

-204.59 

-204.61 

-204.67 

-204.6 

0.500 

-205.96 

-205.96 

-206.1 

1.000 

-206.04 

-206.09 

-206.1 

2.000 

-205.95 

-206.03 

-206.1 

5.000 

-206.35 

-206.10 

-206.13 

-206.34 

-206.2 

10.000 

-206.34 

-206.18 

-206.35 

-206.3 

The  maximum  spread  of  sensitivity  data  for  the  F50  hydrophone  is 
0.4  dB  over  the  frequencies  covered.  Some  of  this  spre  :d  is  caused  by 
the  increased  effect  of  diffraction  of  the  F50  over  that  of  the  longer 
F37 ,  and  some  is  due  to  the  increased  distance  error  because  the  F50 
weighs  less  and  has  a  lighter  cable  than  the  F37,  which  reduces  its 
tendency  to  hang  straight  when  freely  suspended.  The  diffraction  effect 
is  estimated  to  reduce  sensitivity  at  10  kHz  by  from  0.15  to  0.20  dB 
below  that  at  lower  frequencies.  At  500  Hz  the  signal-to-noise  ratio  was 
so  low  that  the  estimates  were  judged  to  be  biased  by  about  +0.1  dB. 

Taking  the  10  estimates  below  10  kHz  as  being  estimates  for  the  same  true 
value,  one  obtains  -206.11  dB  re  1  V/pFa  for  their  average.  The  standard 
deviation  of  these  points  is  0.13  dB,  making  the  standard  error  of  the 
average  to  be  about  0.04  dB.  Tables  for  Student’s  t-distribution  indi¬ 
cate  that  the  confidence  interval  is  -206.21  to  -206.01  dB  with  a  con¬ 
fidence  coefficient  of  0.95.  Again,  because  of  the  possibility  of  posi¬ 
tive  bias,  the  interval  -206.26  to  -206.01  dB  T-dth  95%  certainty  is  not 
unreasonable. 

In  view  of  -ne  planned  design  of  a  new  calibration  system  for  the 
Leesburg  Facility,  it  seems  appropriate  to  say  what  probability  interval 
may  be  expected  from  reciprocity  measurements  based  on  the  methods 
employed  here.  One  would  expect  a  probability  interval  of  about  ±0.30  dB 
to  have  a  probability  of  0.95  of  containing  the  mean  value  of  hydrophone 
sensitivity  for  the  F50  with  comparable  values  of  about  ±0.15  dB  and 
0.95  for  the  F37.  These  intervals  can  be  reduced  by  closer  attention  to 
signal-to-noise  ratios,  environmental  conditions,  and  distance  measure¬ 
ments.  By  using  the  sing-around  velocimeter  principle  to  measure  distance 
at  the  depth  of  the  suspended  transducers  instead  of  at  the  top  of  the 
rigging,  it  may  be  possible  to  achieve  smaller  probability  intervals. 
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A  Preliminary  Study  of  the  Precision  and  Accuracy  of  Hydrophone 
Calibration  at  the  Established  USRD  Facilities 

In  this  section,  the  measurements  of  hydrophone  sensitivity  made  at 
the  Lees-burg  Facility,  the  Anechoic  Tank  Facility,  System  Y.  Tube  Facility, 
Center  Pier  Facility,  Digital  Measuring  System  Facility,  and  North  Pier 
Facility  are  presented.  To  whatever  extent  is  possible,  statistical  state¬ 
ments  about  errors  are  based  on  a  very  limited  amount  of  data  and  on  the 
use  of  the  average  values  from  the  calibrations  of  the  F37  and  the  F50  as 
"true"  values.  Also,  estimates  at  the  different  frequencies  were  assumed 
independent  for  a  single  run,  even  though  the  transducers  were  not  rerigged. 

Lessburg  Calibration  Facility 

Three  calibrations  each  were  made  on  the  F37  serial  A40  and  the  F50 
serial  21  at  the  depth  1480  cm  and  separation  100  cm.  These  were  compari¬ 
son  calibrations  with  the  F37  serial  5  (previously  calibrated  by  reci¬ 
procity)  used  aS  the  standard.  Table  III  shows  the  results  of  these 
calibrations  and  their  comparisons  with  the  "true"  values.  In  this  table, 

B  denotes  the  difference  between  estimated  sensitivity  and  "true" 
sensitivity. 

Table  III.  Free-field  voltage  sensitivity  M  (dB  re  1  V/uPa) 
and  bias  estimate  B  for  two  hydrophones  measured  at  the 
Leesburg  Facility. 


Hydro¬ 

phones 

Freq 

(kHz) 

A 

Mi 

Bl 

m2 

B2 

m3 

h 

F37 

0.5 

-204.3 

+0.4 

-204.5 

+0.2 

-204.5 

+0.2 

1.0 

-204.3 

+0.4 

-204.6 

+0.1 

-204.5 

+0.2 

2.0 

-204.3 

+0.4 

-204.6 

+0.1 

-204.4 

+0.3 

5.0 

-204.3 

+0.4 

-204.5 

+0.2 

-204.4 

+0.3 

10.0 

-204.4 

+0.2 

-204.7 

-0.1 

-204.6 

0 

F50 

0.5 

-205.8 

+0.3 

-205.8 

+0.3 

-205.9 

+0.2 

1.0 

-205.8 

+0.3 

-205.9 

+0.2 

-205.9 

+0.2 

2.0 

-205.7 

+0.4 

-205.9 

+0.2 

-205.9 

+0.2 

5.0 

-205.8 

+0.3 

-206.0 

+0.1 

-206.0 

+0.2 

10.0 

-206.0 

+0.3 

-206.1 

+0.1 

-206.2 

+0.1 

The  maximum  deviation  of  these  results  from  the  "true"  value  is 
+0.4  dB.  The  estimated  bias  B  for  the  F37  serial  A40  is  +0.22  dB.  If  one 
assumes  distance  error  is  negligible,  then  there  are  15  independent  esti¬ 
mates  on  each  hydrophone.  Maximum  deviation  on  repeatability  is  0.3  dB. 

The  estimated  standard  error,  about  B  for  the  F37,  is  o-  =  0.05  dB,  which, 

B 

using  the  Student  t-distribution,  gives  the  confidence  interval  for  the 
bias  as  approximately  +0.12  <  B  <  +0.32  dB  with  95%  certainty.  This 
amount  of  bias  can  be  due  to  a  bias  in  the  calibration  of  the  F37  serial  5 
used  as  a  standard.  If  one  assumes  this  to  be  the  case  and  if  comparison 
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were  made  using  an  unbiased  standard,  confidence  limits  at  the  Leesburg 
Facility  for  calibrating  an  F37  would  be  about  ±0.4  dB  with  95%  certainty, 
assuming  a  normal  distribution  and  that  au  -  a. 

mm 

The  estimated  bias  for  the  F50  serial  21  is  -J-0 .23  dB.  Again  assuming 
negligible  distance  error,  the  estimated  standard  error  about  B  is  0.01  dB 
which,  using  Student's  t-distribution,  gives  the  confidence  interval  for 
the  bias  as  0.21  <  B  <  0.25  dB  with  95%  certainty.  If  the  biases  were  due 
to  calibration  of  the  F37  serial  5  alone,  then  comparison  with  an  unbiased 
standard  would  give  confidence  limits  of  about  ±0.1  with  95%  certainty. 

In  calibrating  both  the  F37  serial  A40  and  the  F50  serial  21,  the  esti¬ 
mated  bias  was  about  0.2  dB.  The  source  of  the  bias  was  not  ascertained, 
but  would  seem  to  be  either  a  distance  bias  or  bias  in  the  sensitivity  of 
the  standard  used  in  the  comparison. 

Anechoic  Tank  Calibration  Facility 

Each  of  the  three  calibrations  made  by  the  Anechoic  Tank  Facility  on 
the  two  test  hydrophones  are' shown  in  Table  IV.  These  measurements  were 
made  by  comparison  with  a  "known"  standard. 

Table  IV.  Free-field  voltage  sensitivity  A  (dB  re  1  V/yPa) 
and  bias  estimate  6  for  two  hydrophones  measured  in  the 
Anechoic  Tank  Facility 


Hydro-  Freq 

phones  (kHz)  Mi  Bj  M2  §2  M3  §3 


F37  0.5 

1.0 
2.0 
5.0 
10.0 

F50  0.5 

1.0 
2.0 
5.0 
10.0 


-204.8 

-0.1 

-204.8 

-0.1 

-204.8 

-0.1 

-204.8 

-0.1 

-204.8 

-0.2 

-206.2 

-0.1 

-206.1 

0 

-206.1 

0 

-206.2 

0 

-206.2 

+0.1 

-204.5 

+0.2 

-204.5 

+0.2 

-204.5 

+0.2 

-204.5 

+0.2 

-204.5 

+0.3 

-206.0 

+0.1 

-206.0 

+0.1 

-206.0 

+0.1 

-206.0 

+0.2 

-206.0 

+0.3 

-204.5 

+0.2 

-204.5 

+0.2 

-204.5 

+0.2 

-204.5 

+0.2 

-204.5 

+0.3 

-206.0 

+0.1 

-206.0 

+0.1 

-206.0 

+0.1 

-206.0 

+0.2 

-206.0 

+0.3 

The  maximum  deviation  from  "true"  value  on  these  measurements  is 
+0.3  dB.  Estimated  biases  for  the  F37  and  F50  are  +0.08  and  +0.11  dB, 
respectively.  Maximum  deviation  from  repeatability  is  0.5  dB,  while  the 
estimated  standard  errors  about  the  B's  are  0.05  dB  for  the  F37  and 
0.03  dB  for  the  F50.  These  values  correspond  to  confidence  intervals  of 
-0.03  <  B  <  +0.19  dB  and  +0.05  <  B  <  +0.17  dB,  respectively,  at  the  95% 
confidence  level.  The  F37  cannot  be  said  to  have  a  bias  in  its  estimated 
sensitivity,  but  the  F50  can  at  the  95%  level  of  certainty.  Confidence 
intervals  on  estimated  hydrophone  sensitivities  calibrated  in  the  Anechoic 


23 


V 


Tank  are  about  M  -  0.1  <  M  <  M  +  0.3  dB  for  the  F37  and  M  +  0.0  <  M  < 

M  +  0.2  dB  for  the  F50  to  about  95%  certainty. 

Low  Frequency  Facility,  System  K 

Measurements  made  in  System  K  on  the  two  standard  hydrophones  were 
limited  to  frequencies  at  or  below  1  kHz.  Comparison  is  made  to  refer¬ 
ence  values  of  -204.7  dB  re  1  V/  Pa  for  the  F37  and  -206.1  dB  re  1  V/yPa 
for  the  F50.  Results  of  these  measurements  are  shown  in  Table  V.  These 
measurements  were  made  by  comparison  with  a  small  hard  hydrophone  that 
previously  had  been  calibrated  in  a  coupler. 


Table  V.  Free-field  voltage  sensitivity  M  (dB  re  lyV/  Pa) 
and  bias  estimate  B  for  two  hydrophones  measured  in  the 
Low  Frequency  Facility,  System  K. 


Hydro¬ 

phones 

Freq 

(Hz) 

A 

Ml 

Bl 

m2 

A 

b2 

m3 

A 

B3 

F37 

10 

-204.7 

0 

-204.8 

-0.1 

-204.7 

0 

20 

-204.7 

0 

-204.8 

-0.1 

-204.7 

0 

50 

-204.7 

0 

-204.8 

-0.1 

-204.7 

0 

100 

-204.7 

0 

-204.8 

-0.1 

-204.7 

0 

200 

-204.7 

0 

-204.8 

-0.1 

-204.7 

0 

500 

-204.6 

+0.1 

-204.8 

-0.1 

-204.6 

+0.1 

800 

-204.6 

+0.1 

-204.8 

-0.1 

-204.6 

+0.1 

1000 

-204.6 

+0.1 

-204.8 

-0.1 

-204.6 

+0.1 

F50 

10 

-205.8 

+0.3 

-205.8 

+0.3 

-205.8 

+0.3 

20 

-205.8 

+0.3 

-205.8 

+0.3 

-205.8 

+0.3 

50 

-205.8 

+0.3 

-205.8 

+0.3 

-205.8 

+0.3 

100 

-205.8 

+0.3 

-205.8 

+0.3 

-205.9 

+0.2 

200 

-205.9 

+0.2 

-205.8 

+0.3 

-205.9 

+0.2 

500 

-205.8 

+0.3 

-205.8 

+0.3 

-206.0 

+0.1 

800 

-205.9 

+0.2 

-205.8 

+0.3 

-206.1 

0 

1000 

-205.9 

+0.2 

-205.8 

+0.3 

-205.2 

-0.1 

The  maximum  deviations  from  "true"  values  on  the  F37  and  F50  are 
±0.1  dB  and  +0.3  dB,  respectively.  Estimated  bias  B  is  -0.01  dB  for  the 
F37  and  +0.24  for  the  F50.  Maximum  deviation  from  repeatability  is 
0.3  dB,  while  the  estimated  standard  errors  about  the  B's.are  0.03  dB  for 
the  F37  and  0.04  for  the  F50.  These  values  correspond  to  confidence 
intervals  of  -0.02  <  B  <  0.00  for  the  F37  and  +0.22  <  B  <  +0.26  for  the 
F50  with  95%  certainty.  Confidence  intervals  on  the  estimated  hydrophone 
sensitivities  calibrated  with  System  K  are  about  M  -  0.1  <M<$+0.1dB 
for  the  F37  and  about  M+0.2  <M  <  M  +  0.3  dB  for  the  F50  with  about  95% 
certainty . 
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Lake  Gem  Mar y  Facility,  Center  Pier  System 

Results  of  measurements  made  by  comparison  for  the  two  standard  hydro¬ 
phones  are  shown  in  Table  VI. 


Table  VI.  Free-field  voltage  sensitivity  M  (dB  re  1  V/yPa) 
and  bias  estimate  B  for  two  hydrophones  measured  in  the 
Lake  Gem  Mary  Facility,  Center  Pier  System. 


Hydro¬ 

phone 

Freq 

(kHz) 

Mi 

A 

Bl 

m2 

A 

B2 

m3 

A 

B3 

F37 

0.5 

-205.0 

-0.3 

-204.5 

+0.2 

-205.2 

-0.5 

1.0 

-205.0 

-0.3 

-204.5 

+0.2 

-205.2 

-0.5 

2.0 

-204.9 

-0.2 

-204.6 

+0.1 

-204.9 

-0.2 

5.0 

-205.0 

-0.3 

-204.5 

+0.2 

-205.0 

-0.3 

10.0 

-204.6 

0 

-204.5 

+0.1 

-204.8 

-0.2 

F50 

0.5 

-206.0 

+0.1 

-206.0 

+0.1 

-206.4 

-0.3 

1.0 

-206.0 

+0.1 

-206.0 

+0.1 

-206.4 

-0.3 

2.0 

-206.1 

0 

-206.0 

+0.1 

-2C6.3 

-0.2 

5.0 

-206.3 

-0.1 

-206.3 

-0.1 

-206.0 

+0.2 

10.0 

-206.2 

+0.1 

-206.3 

0 

-206.2 

+0.1 

The  maximum  deviations  from  "true"  values  on  the  F37  and  F50  are 
-0.5  dB  and  -0.3  dB,  respectively.  Estimated  bias  B  is  -0.13  for  the  F37 
and  -0.07  for  the  F50.  Maximum  deviation  from_repeatability  is  0.7  dB, 
while  the  estimated  standard  errors  about  the  B's  are  0.06  dB  for  the  F37 
and  0.05  dB  for  the  F50.  These  values  correspond  to  confidence  intervals 
of  -0.25  <  B  <  -0.01  dB  and  -0.18  <  B  <  +0.04  dB  for  the  F37  and  F50, 
r aspectively ,  with  95%  certainty.  Confidence  intervals  on  the  estimated 
hydrophone  sensitivities  calibrated  with  the  center  pier  system  are  about 
M  -  0.4  <  M  <  A  +  0.1  dB  for  the  F37  and  about  K  -  0.3  <  M  <  M  +  0.1  dB 
for  the  F50  with  about  95%  certainty. 

Lake  Gem  Mary  Facility,  Digital  Measuring  System 

Results  of  measurements  made  by  comparison  for  the  two  standard  hydro¬ 
phones  are  shown  in  Table  VII. 

The  maximum  deviation  from  "true"  values  on  the  F37  and  F50  are 
-1.0  dB  and  -0.6  dB  respectively.  Estimated  bias  B  is  -0.72  dB  and 
-0.19  dB  for  the  F37  and  F50,  respectively.  Maximum  deviation  from 
repeatability  is  0.6  dB,  while  the  estimated  standard  errors  about  the  B's 
are  0.03  dB  for  the  F37  and  0.06  dB  for  the  F50.  These  values  correspond 
to  confidence  intervals  of  -0.78  <  B  <  -0.66  and  -0.31  <  B  <  -0.07  dB  for 
the  F37  and  F50,  respectively,  with  95%  certainty.  Confidence  intervals 
on  the  estimated  hydrophone  sensitivities  calibrated  with  the  digital 
measuring  system  are  about  M  -  1.0  <  M  <  M  -  0.5  dB  for  the  F37  and  about 
A  -  0.7  <M<A+0.2dB  for  the  F50  with  about  95%  certainty. 
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Table  VII.  Free-field  voltage  sensitivity  M  (dB  re  3  V/yPa) 
and  bias  estimate  B  for  two  hydrophones  measured  in  the 
Lake  Gem  Mary  Facility,  Digital  Measuring  System. 


Hydro¬ 

phone 

Freq 

(kHz) 

A 

Mi 

A 

Bl 

CM 

<s 

A 

B2 

M3 

A 

B3 

F37 

0.5 

-205.4 

-0.7 

-205.3 

-0.6 

-205.5 

-0.8 

1.0 

-205.5 

-0.8 

-205.5 

-0.8 

-205.5 

-0.8 

2.0 

-205.4 

-0.7 

-205.3 

-0.6 

-205.7 

-1.0 

5.0 

-205.5 

-0.8 

-205.3 

-0.6 

-204.7 

-0.6 

10.0 

-205.2 

-0.6 

-205.4 

-0.8 

-204.6 

-0.6 

F50 

0.5 

-206.3 

-0.2 

-206.4 

-0.3 

-206.3 

-0.2 

1.0 

-206.1 

-0.1 

-206.5 

-0.4 

-206.3 

-0.2 

2.0 

-206.2 

-0.1 

-206.4 

-0.3 

-206.3 

-0.2 

5.0 

-206.4 

-0.2 

-206.3 

-0.1 

-205.8 

+0.4 

10.0 

-206.7 

-0.4 

-206.9 

-0.6 

-206.3 

0 

Lake  Gem  Mary  Facility ,  North  Pier  System 

Measurements  made  on  the  North  Pier  System  showed  that  the  system  was 
not  suitable  for  calibrating  at  that  time.  No  calibrations  were  being  per¬ 
formed  on  the  system  at  that  time  because  facility  personnel  felt  the 
system  was  inaccurate.  Measurements  showed  a  -1.0  to  -1.5  dB  bias  in  cali¬ 
brating  the  F37  serial  A40  and  F50  serial  21  with  some  values  as  high  as 
-1.8  dB  off.  The  system  has  been  repaired  and  is  said  to  be  working 
properly,  but  no  new  measurements  have  been  made  on  it  for  this  report. 


SUMMARY  AND  DISCUSSION  OF  SYSTEM  ERROR  MEASUREMENTS 

Figure  9  is  a  precision  and  accuracy  diagram  for  all  of  the  systems 
for  both  the  F37  serial  A40  and  the  F50  serial  21  measurements  relative  to 
the  estimates  obtained  with  the  simple  calibration  system  as  the  "true" 
value.  The  estimated  bias  |b|  =  |m  -  "M" |  is  shown  for  each  hydrophone  as 
calculated  from  measurements  made  on  each  system.  Also  shown  are  the 
confidence  intervals  associated  with  the  "true"  value  "M" .  The  upper  of 
the  pair  of  confidence  intervals  represents  the  95%  certainty  for  M.  The 
lower  confidence  interval  is  for  M  of  the  probability  that  95%  of  the 
estimates  for  hydrophone  sensitivity  made  from  single  measurements  will 
fall  within  the  marked  interval,  assuming  aM  =  aM  and  a  normal  distribu¬ 
tion  for  M.  The  X's  indicate  values  that  fell  outside  of  the  confidence 
interval,  while  the  O's  indicate  other  extreme  departures  from  M.  Out  of 
150  values,  only  6  fell  outside  the  confidence  intervals  (about  96%  were 
within  the  intervals) . 

From  these  measurements,  one  can  conjecture  that  the  USRD  is  very  well 
justified  in  claiming  that  an  interval  of  ±1  dB  about  the  true  value  of 
hydrophone  sensitivity  will  contain  at  least  95%  of  the  estimates  of 
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Fig.  9.  Precision  and 
accuracy  diagram  for 
hydrophone  sensitivity 
measurements  at  USRD 
facilities:  (a)  Leesburg, 

(b)  Anechoic  Tank, 

(c)  System  K,  (d)  Center 
Pier,  (e)  Digital 
Measuring  System. 


hydrophone  sensitivity  for  small  standard  hydrophones  operated  far  below 
resonance.  At  90%  certainty  all  systems  except  the  digital  measuring 
system  seem  to  have  an  interval  of  ±0.5  dB  associated  with  them.  Care- 
full  attention  to  sources  of  bias,  repetitive  calibration,  and  modern 
signal  processing  techniques  probably  could  lead  to  a  confidence  interval 
of  ±0.1  to  ±0.2  dB  without  much  added  system  complexity. 


RECOMMENDATIONS 

A.  Study  of  system  accuracy  should  continue  with  a  system-by-system 
approach  instead  of  using  an  over-all  approach  in  order  to  obtain  enough 
data  from  a  system  to  be  able  to  pinpoint  random  and  systematic  errors. 

B.  A  study  of  distance  errors  associated  with  rigging  of  transducers 
should  be  undertaken.  This  could  be  done  with  sufficient  accuracy  by- 
application  of  the  sing-around  velocimeter  principle. 

C.  A  study  of  signal  processing  techniques  aimed  at  calibration  with 
low  signal-to-noise  ratios  is  recommended.  This  could  lead  to  better 
free-field  calibration  at  low  frequencies  where  noise  is  high  and  where 
difficulty  is  encountered  in  building  small,  high-intensity  sound  sources. 

D.  A  study  of  the  radiation  impedance  of  standard  hydrophones  is 
recommended  with  a  view  toward  developing  cutoff  filters  for  some  stand¬ 
ards,  which  would  limit  their  use  to  the  frequency  range  where  they  can 
be  calibrated  and  used  without  too  much  attention  to  boundary  conditions. 
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APPENDIX  A 


This  appendix  consists  of  photographs,  drawings,  and  typical  sensi¬ 
tivity  curves  for  the  F37  and  F50  hydrophones. 


Fig.  Al.  F37  hydrophone 
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Fig.  A2.  Nominal  dimensions  of  F37 
hydrophone . 


Pig.  A3.  F50  hydrophone 
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Pig.  A4.  Nominal  dimensions  of  F50 
hydrophone . 


